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ABSTRACT 



RARE KAON DECAYS AND CP VIOLATION 
SEPTEMBER 1997 
FABRIZIO GABBIANI 
LAUREA IN FISICA, UNIVERSITA DEGLI STUDI DI PADOVA 
Ph.D., UNIVERSITY OF MASSACHUSETTS AMHERST 
Directed by: Professor John F. Donoghue 

Rare kaon decays are an important testing ground of the electroweak flavor 
theory. They can provide new signals of CP-violating phenomena and open a 
window into physics beyond the Standard Model. The interplay of long-distance 
QCD effects in strangeness-changing transitions can be analyzed with Chiral Per- 
turbation Theory techniques. Some theoretical predictions obtained within this 
framework for radiative kaon decays are reviewed, together with the present ex- 
perimental status. In particular, two rare kaon decays are analyzed: The first 
decay, vr^e+e", is being searched for as a signal of direct AS" — 1 CP vio- 

lation. We provide a thorough updating of the analysis of the three components 
of the decay: 1) Direct CP violation, 2) CP violation through the mass matrix 
and 3) CP-conserving (two-photon) contributions. First the chiral calculation of 
the Ks — > 7r°e+e~ rate, due to Ecker, Pich and de Rafael, is updated to include 

vii 



recent results on the nonleptonic amplitude. Then we systematically explore the 
uncertainties in this method. These appear to be so large that they will obscure 
the direct CP violation unless it is possible to measure the Kg — > 7r°e'^e~ rate. The 
CP-conserving amplitude remains somewhat uncertain, but present indications are 
such that there may be a sizable CP- violating asymmetry in the e"'",e~ energies 
from the interference of CP-conserving and CP- violating amplitudes and this may 
potentially be useful in determining whether direct CP violation is present. The 
second decay, Kl 7r^je'^e~, which occurs at a higher rate than the nonradiative 
process Kl — > 7r°e+e~, can be a background to CP violation studies using the latter 
reaction. It also has interest in its own right in the context of Chiral Perturbation 
Theory, through its relation to the decay Kl — > 77^77. The leading order chiral 
loop contribution to Kl — > 7r°7e'*"e~, including the (ge+ + ?e-)^/^^ dependence, 
is completely calculable. We present this result and also include the higher order 
modifications that are required in the analysis of Kl — > 77*^77. 
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CHAPTER 1 



INTRODUCTION 



High-precision experiments on rare kaon decays offer tlie possibiiity of unrav- 
eling new pfiysics beyond tfie Standard Model. Searching for forbidden flavor- 
changing processes [1] at the 10"^° level [BR{Kl /xe) < 3.3 x 10-^^ BR{Kl 
vrVe) < 3.2 x 10"^ BR{K+ 7r+/xe) < 2.1 x 10-^°, ...], one is actually ex- 
ploring energy scales above the 10 TeV region. The study of allowed (but highly 
suppressed) decay modes provides, at the same time, very interesting tests of the 
Standard Model itself. Electromagnetic-induced nonleptonic weak transitions and 
higher-order weak processes are a useful tool to improve our understanding of the 
interplay among electromagnetic, weak and strong interactions. In addition, new 
signals of CP violation, which would help to elucidate the source of CP-violating 
phenomena, can be looked for. 

In this chapter we shall briefly describe a few relevant rare kaon decays, reserving 
the rest of this thesis for the processes Kl — > 7r°77, Kl — > 7r°e+e~ and Kl 
7r°7e+e", together with the related decays tt'^jj and — > ir'^e'^e'. 

Since the kaon mass is a very low energy scale, the theoretical analysis of nonlep- 
tonic kaon decays is highly non-trivial. While the underlying flavor-changing weak 
transitions among the constituent quarks are associated with the W-mass scale. 
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the corresponding hadronic amplitudes are governed by the long-distance behavior 
of the strong interactions, i.e. the confinement regime of QCD. 

The standard short-distance approach to weak transitions makes use of the 
asymptotic freedom property of QCD to successively integrate out the fields with 
heavy masses down to scales /i < rric- Using the operator product expansion (OPE) 
and renormalization-group techniques, one gets an effective AS — 1 hamiltonian, 

^f/=' = EQiM + h.c., (1.1) 

which is a sum of local four-fermion operators Qi, constructed with the light de- 
grees of freedom {u, d, s; e, ^u, ui), modulated by Wilson coefficients Ci(/i) which are 
functions of the heavy {W, t, b, c, r) masses. The overall renormalization scale /i 
separates the short- (M > n) and long- (m < n) distance contributions, which 
are contained in Ci(/x) and Qi, respectively. The physical amplitudes are of course 
independent of thus, the exphcit scale (and scheme) dependence of the Wilson 
coefficients, should cancel exactly with the corresponding dependence of the Qi 
matrix elements between on-shell states. 

Our knowledge of the AS — 1 effective hamiltonian has improved considerably 
in recent years, thanks to the completion of the next-to-leading logarithmic order 
calculation of the Wilson coefficients [2]. All gluonic corrections of 0{a^t^) and 
0{a'^'^H'^) are already known, where t = log (M/m) refers to the logarithm of any 
ratio of heavy-mass scales {M,m > ii). Moreover, the full rrit/Mw dependence (at 
lowest order in ag) has been taken into account. 

Unfortunately, in order to predict the physical amplitudes one is still confronted 
with the calculation of the hadronic matrix elements of the quark operators. This 
is a very difficult problem, which so far remains unsolved. We have only been 
able to obtain rough estimates using different approximations (vacuum saturation. 
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Nc — > oo limit, QCD low-energy effective action, . . . ) or applying QCD techniques 
(lattice, QCD sum rules) which suffer from their own technical limitations. 

Below the resonance region (// < mp) the strong interaction dynamics can be 
better understood with global symmetry considerations. We can take advantage of 
the fact that the pseudoscalar mesons are the lowest energy modes of the hadronic 
spectrum: They correspond to the octet of Goldstone bosons associated with the 
dynamical chiral symmetry breaking of QCD, SU{?>)l ® SU{3)r SU{3)v- The 
low-energy implications of the QCD symmetries can then be worked out through 
an effective lagrangian containing only the Goldstone modes. The effective Chiral 
Perturbation Theory [3] (ChPTh) formulation of the Standard Model is an ideal 
framework to describe kaon decays [4] . This is because in K decays the only physical 
states that appear are pseudoscalar mesons, photons and leptons, and because the 
characteristic momenta involved are small compared to the natural scale of chiral 
symmetry breaking (A^^ ~ 1 GeV). 

1.1 K ^ TTiyV 

The decays K — > ttw proceed through flavor-changing neutral current effects. 
These arise in the Standard Model only at second (one-loop) order in the elec- 
troweak interaction (Z-penguin and W-box diagrams. Fig. 1) and are additionally 
GIM suppressed. 

The branching fractions are thus very small, at the level of 10~^°, which makes 
these modes rather challenging to detect. However, K — > nuV have long been 
known to be reliably calculable, in contrast to most other decay modes of interest. 
A measurement of K'^ — > tt'^uT' and — > tt^i'T' will therefore be an extremely 
useful test of flavor physics. Over the recent years important reflnements have 



a) 



b) 
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Figure 1: Diagrams inducing the decay s — > dvv in the Standard Model. 

been added to the theoretical treatment oi K ^ tivV. Let us briefly summarize 
the main aspects of why K ttuV is theoretically so favorable and what recent 
developments have contributed to emphasize this point. 

a) First, K — > TruV is semileptonic. The relevant hadronic matrix elements, such 

as {7i\{sd,)v-A\K) , are just matrix elements of a current operator between ha- 
dronic states, which are already considerably simpler objects than the matrix 
elements of four-quark operators encountered in many other observables {K — 
K mixing, e'/e). Moreover, they are related to the matrix element 



by isospin symmetry. The latter quantity can be extracted from the well 



{A{su)v.a\K+) 



(1.2) 



measured leading semileptonic decay K'^ — > tt'^Iv. Although isospin is a fairly 
good symmetry, it is still broken by the small up-down quark mass difference 
and by electromagnetic effects. These manifest themselves in differences of the 
neutral versus charged kaon (pion) masses (affecting phase space), corrections 
to the isospin limit in the form factors and electromagnetic radiative effects. 
Marciano and Parsa [5] have analyzed these corrections and found an overall 
reduction in the branching ratio by 10% for n^uV and 5.6% for Kl — > 

b) Long distance contributions are systematically suppressed as 0{AQQ]^/ml) com- 

pared to the charm contribution (which is part of the short distance ampli- 
tude). This feature is related to the hard (~ m^) GIM suppression pattern 
shown by the Z-penguin and W-box diagrams, and the absence of virtual 
photon amplitudes. Long distance contributions have been examined quanti- 
tatively [6] and shown to be numerically negligible (below ~ 5% of the charm 
amplitude) . 

c) The preceding discussion implies that K tti/V are short distance dominated by 

top- and charm-loops in general. The relevant short distance QCD effects can 
be treated in perturbation theory and have been calculated at next-to-leading 
order [7]. This allowed to substantially reduce (for K'^) or even practically 
eliminate (-R'l) the leading order scale ambiguities, which are the dominant 
uncertainties in the leading order result. 

The decay amplitude for K — > ttuu, 

i=c,t 

(1-3) 
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involves the hadronic matrix element of the AjS" = 1 vector current, which (assuming 
isospin symmetry) can be obtained from Ki^ decays. In the ChPTh framework, the 
needed hadronic matrix element is known at this allows one to estimate the 

relevant isospin- violating corrections reliably [8, 5]. 

Summing over the three neutrino flavors and expressing the quark-mixing fac- 
tors through the Wolfenstein parameters [9] A, A, p and 77, one can write the 
approximate formula [2]: 



BR{K+ TT+uV) ^ 1.93 X \r]^ + {po - p) 



p^^lA. (1.4) 



The departure of po from unity measures the impact of the charm contribution. 

With the presently favored values for the quark-mixing parameters, the branch- 
ing ratio is predicted to be in the range [2] 



BR{K+ n+uV) = (9.1 ± 3.2) x 10"^^ , (1.5) 

to be compared with the present experimental upper bound [10] BR{K'^ — > t:'^vV) 
< 2.4x10-9 (90% CL). 

What is actually measured is the transition tt"*"-!- nothing; therefore, 

the experimental search for this process can also be used to set limits on possible 
exotic decay modes like 7r+X°, where denotes an undetected light Higgs 

or Goldstone boson (axion, familon, majoron, . . .). 

The CP- violating decay Kl — > t^^vV has been suggested [11] as a good candidate 
to look for pure direct CP-violating transitions. The contribution coming from 
indirect (mass- matrix) CP violation via mixing is very small [11]: BR{Kl — > 

t^^vV)mm ~ 5 xlO"^^. The decay proceeds almost entirely through direct CP 
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Table 1: Important properties and results for K — > nvv. 











CP-conserving 


CP-violating 


CKM 


Vu 


~ Jcp ~ ?7 


scale uncert. (BR) 


±20% (LO) ^ ±5% (NLO) 


±10% (LO) ^ ±1% (NLO) 


BR (SM) 


(0.9 ± 0.3) X R)-^" 


(2.8 ± 1.7) X IQ-^^ 


exp. limit 


< 2.4 X 10-9 BNL 787 [10] 


< 5.8 X 10"^ FNAL 799 [12] 



violation, and is completely dominated by short-distance loop diagrams with top 

quark exchanges [2]: 

BR(Kl 7r%v) ^ 8.07 x 10"^^ 77^ x]-^^ . (1.6) 

The present experimental upper bound [12], BR{Kl tt^uV) < 5.8 x 10"^ (90% 
CL) , is still far away from the expected range [2] 

BR{Kl = (2.8 ± 1.7) x 10"^^ . (1.7) 

Nevertheless, the experimental prospects to reach the required sensitivity in the 
near future look rather promising. The clean observation of just a single unam- 
biguous event would indicate the existence of CP-violating A^" = 1 transitions. 

In Table 1 we have summarized some of the main features of — > 7r~^i>V and 

While already K'^ — > n'^uV can be reliably calculated, the situation is even 
better for Kl t^^i^V. The charm sector, in ti^vV the dominant source of 

uncertainty, is completely negligible for Kl — > t^^vV (0.1% effect on the branching 
ratio). Long distance contributions ( ^ 0.1%) and also the indirect CP violation 
effect ( ^ 1%) are likewise negligible. In summary, the total theoretical uncer- 
tainties, from perturbation theory in the top sector and in the isospin breaking 
corrections, are safely below 2 — 3% for BR{Kl — > t^^vV). This makes this decay 
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mode truly unique and very promising for phenomenological applications. Note 
that the range given as the Standard Model prediction in Table 1 arises from our, 
at present, limited knowledge of Standard Model parameters (CKM), and not from 
intrinsic uncertainties in calculating BR[Kl ti^uV). 

With a measurement of BR{K~^ tt'^i/V) and BR{Kl tt^vV) available very 
interesting phenomenological studies could be performed. For instance, BR{K^ — > 
Ti^uV) and BR{Kl n^vV) together determine the unitarity triangle (Wolfenstein 
parameters p and rf) completely. 

The expected accuracy with ±10% branching ratio measurements is comparable 
to the one that can be achieved by CP violation studies at B factories before the 
LHC era [13]. The quantity BR{Kl n^uv) by itself offers probably the best 
precision in determining ImV^*T4(i or, equivalently, the Jarlskog parameter [14] 



The charged mode K'^ 'k'^vV is being currently pursued by Brookhaven 
experiment E787. The latest published result [10] gives an upper limit which is 
about a factor 25 above the Standard Model range. Several improvements have 
been implemented since then and the SM sensitivity is expected to be reached 
in the near future [15]. Recently an experiment has been proposed to measure 
— > 'K^vV at the Fermilab Main Injector [16]. Concerning — > 'k^vV^ a 
proposal exists at Brookhaven (BNL E926) to measure this decay at the AGS with 
a sensitivity of 0(10"^^) (see [15]). There are furthermore plans to pursue this 
mode with comparable sensitivity at Fermilab [17] and KEK [18]. 




(1.8) 





Figure 2: Feynman diagrams for — > 77. 
1.2 Kg — ^ 77 

The symmetry constraints do not allow any direct tree-level i^i77 coupling 
at 0{E'^) {K^2 ^^^^^ CP-even and CP-odd eigenstates, respectively). This 

decay proceeds then through a loop of charged pions and kaons as shown in Fig. 2. 
Since there are no possible counter-terms to renormalize divergences, the one-loop 
amplitude is necessarily finite. Although each of the four diagrams in Fig. 2 is 
quadratically divergent, these divergences cancel in the sum. The resulting pre- 
diction [19], BR{Ks — > 77) = 2.0 X 10~^, is in very good agreement with the 
experimental measurement [20]: 
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BR{Ks 77) = (2.4 ± 0.9) x 10 



-6 



(1.9) 



1.3 Kl,s fJ^^I^ 



There are well-known short-distance contributions [2] (electroweak penguins 
and box diagrams) to the decay A*^A*~- This part of the amplitude is 

sensitive to the Wolfenstein parameter p. However, this transition is dominated 
by long-distance physics. The main contribution proceeds through a two-photon 
intermediate state: K2 — > 7*7* — > jjf^jJT . Contrary to — > 77, the prediction 
for the K2 "PI decay is very uncertain, because the first non-zero contribution 
occurs^ at 0{E^^. That makes very difficult any attempt to predict the A*'*'/^"' 
amplitude. 

The long distance amplitude consists of a dispersive (Adis) and an absorptive 
contribution (>labs)- The branching fraction can thus be written 



associated with the 77 intermediate state, and knowing BR{Kl — > 77) it is possible 

to extract l^absP = (6.8 ± 0.3) x 10~^ [21]. A^^s on the other hand cannot be 

calculated accurately at present and the estimates are strongly model dependent 

^At 0{E'^), this decay proceeds through a tree-level K2 7r°, rj transition, followed by 7r°, r] 
77 vertices. Because of the Gell-Mann-Okubo relation, the sum of the tt" and r] contributions 
cancels exactly to lowest order. The decay amplitude is then very sensitive to SU{3) breaking. 



BR{Kl ii-^ii-) = \AsB + AdisI' + |Aabs 



(1.10) 



Using the well-known unitarity bound 




(1.11) 



11 




II- 



Figure 3: Feynman diagram for the decay . The K^^*'y* 

vertex is generated through the one-loop diagrams shown 
in Fig. 2. 

[22]. This is rather unfortunate, in particular since BR(Kl — > jJL^jJL~), unlike most 
other rare decays, has already been measured, and this with very good precision: 



BR{Kl ^ = { 



(1.12) 



(6.9 ± 0.4) X 10-9, BNL 791 [23], 
(7.9 ± 0.7) X 10"^ KEK 137 [24]. 

For comparison we note that BR{Ki Ai^At-)sD = (1-3 ± 0.6) x 10"^ is the 
expected branching ratio in the Standard Model based on the short-distance con- 
tribution alone. Due to the fact that A^^^ is largely unknown, — > A*"*"/^" is 
at present not a very useful constraint on CKM parameters. Some improvement 
might be expected from measuring the decay ix'^ix~e^e~ , which could lead to 

a better understanding of the Kl — > 7*7* vertex. First results obtained at Fermilab 
(E799) give BR{Kl jj+u'e+e-) = (2.9l|I) x 10"^- 

The situation is completely different for the Kg decay. A straightforward chiral 
analysis [25] shows that, at lowest order in momenta, the only allowed tree-level 
K^^^li~ coupling corresponds to the CP-odd state Therefore, the 

li'^li" transition can only be generated by a finite non-local two-loop contribution, 
illustrated in Fig. 3. 



12 



The explicit calculation [25] gives: 



3^^^ ^ 2 X 10-, ^i'fr'V - 8 X (1-13) 

well below the present (90% CL) experimental upper limits [26]: BR{Ks —>■ 
< 3.2 X 10-^ BR{Ks e+e-) < 2.8 x 10"^. Although, in view of the smallness 
of the predicted ratios, this calculation seems quite academic, it has important 
implications for CP-violation studies. 

The longitudinal muon polarization Vl in the decay Kl — > /^"'"A*" is an interest- 
ing measure of CP violation. As for every CP-violating observable in the neutral 
kaon system, there are in general two different kinds of contributions to Vl'- Indi- 
rect mass-matrix CP violation through the small admixture of the Kl (e effect), 
and direct CP violation in the K2 — > A*^A*~ decay amplitude. 

In the Standard Model, the direct CP-violating amplitude is induced by Higgs 
exchange with an effective one-loop flavor-changing sdH coupling [27]. The present 
lower bound on the Higgs mass, Mh > 66 GeV (95% CL), implies a conservative 
upper limit I'Pi.dirl < 10~^. Much larger values, Vl ~ O(10~^), appear quite 
naturally in various extensions of the Standard Model [28] . It is worth emphasizing 
that Vl is especially sensitive to the presence of light scalars with CP-violating 
Yukawa couplings. Thus, Vl seems to be a good signature to look for new physics 
beyond the Standard Model; for this to be the case, however, it is very important 
to have a good quantitative understanding of the Standard Model prediction to 
allow us to infer, from a measurement of Vl, the existence of a new CP- violation 
mechanism. 
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The chiral calculation of the — > amplitude allows us to make a reliable 

estimate of the contribution to due to K^-TC^ mixing [25] : 

2 X 10-6 1 



1.9 < \Vl,: 



MM 



X 10'' 



< 2.5. 



;i.l4) 



Taking into account the present experimental errors in BR[Ks 77) and the 
inherent theoretical uncertainties due to uncalculated higher-order corrections, one 
can conclude that experimental indications for \Vl\ > 5 x 10"^ would constitute 
clear evidence for additional mechanisms of CP violation beyond the Standard 
Model. 



1.4 K+ ttVV" 

The rare decay — > 7r+yLi+/x~ has recently been measured at Brookhaven 
(BNL 787) with a branching ratio [29] 

BR{K+ 7r+//+/i~) = (5.0 ± 0.4 ± 0.6) x 10"^ (1.15) 

This compares well with the estimate from ChPTh BR{K+ 77+//+//") = (6.2tg|) 
x 10^^ [30]. The branching ratio is completely determined by the long-distance 
contribution arising from the one-photon exchange amplitude. A short-distance 
amplitude from Z-penguin and W-box diagrams also exists, but is smaller than the 
long-distance part by three orders of magnitude and does therefore not play any role 
in the total rate. However, while the muon pair couples via a vector current (/!//) y 
in the photon amplitude, the electroweak short-distance mechanism also contains 
an axial vector piece {]Ii^)a- The interference term between these contributions is 
odd under parity and gives rise to a parity- violating longitudinal polarization, 
which can be observed as an asymmetry A^^ = (T r — T l) / {T ji + T l) [31, 32]. Tjk^^-) 
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denotes the rate of producing a right- (left-) handed /j,^ in — > 7r+//+//~ decay. 
The effect occurs for a /x~ instead of /x"*" as well, but the polarization measurement 
is much harder in this case. 

Alr is sensitive to the Wolfenstein parameter p. It is a cleaner observable than 
Kl A*"*"/^", although some contamination through long-distance contributions 
cannot be excluded [32]. At any rate, Alr will be an interesting observable to 
study if a sensitive polarization measurement becomes feasible. The Standard 
Model expectation is typically around Alr ~ 0.5%. 

1.5 Summary 

The field of rare kaon decays raises a broad range of interesting topics, ranging 
from ChPTh, over Standard-Model flavordynamics to exotic phenomena, thereby 
covering scales from Aqcd to the weak scale (M^, rrit) and beyond to maybe several 
hundred TeV's. In the present chapter we have focused on the flavor physics of the 
Standard Model and those processes that can be used to test it. Several promising 
examples of short-distance sensitive decay modes exist, whose experimental study 
will provide important clues on flavordynamics. On the theoretical side, progress 
has been achieved over recent years in the calculation of effective Hamiltonians, 
which by now include the complete NLO QCD effects in essentially all cases of 
practical interest. 

The decay Kl — > A*"*"/^" is already measured quite accurately; unfortunately a 
quantitative use of this result for the determination of CKM parameters is strongly 
limited by large hadronic uncertainties. 

The situation looks very bright for K — > nuV. The charged mode is exper- 
imentally already close and its very clean status promises useful results on Vtd- 
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Finally the decay — > 'K^yV is a particular highlight in this field. It could serve 
for instance as the ideal measure of the Jarlskog parameter Jcp. Measuring the 
branching ratio is a real experimental challenge. 

It is to be expected that rare kaon decay phenomena will in the future continue 
to contribute substantially to our understanding of the fundamental interactions. 



CHAPTER 2 



THEORETICAL FRAMEWORK 

2.1 Motivation 

There are three rare decay modes of the long-hved kaon which have interrelated 
theoretical issues: — > 7r'^77, — > 7r'^e+e~ and Kl — > 7r°7e+e~. The first two 
have been extensively studied; the latter has not been previously calculated. It is 
the purpose of this thesis to provide a calculation of the latter two processes and 
describe how they are related to the phenomenology of the first one. 

There is a curious and important inverted hierarchy of these decay modes. 
The rate for the radiative decay Kl — > 7r°7e"'"e~ is a power of a larger than the 
nonradiative transition Kl — > "K^e^e" . This is because the Kl — > 7r°e'^e~ transition 
occurs only through a two-photon intermediate state, or alternatively through a 
one-photon exchange combined with CP violation (which numerically appears to 
be roughly of the same size as the two-photon contribution) [33] . The Kl 7r°e+e~ 
rate is then of order a'^. However, in Kl —>■ 7r°7e+e~ we need only a one-photon 
exchange to the e'^e~, leading to a rate of order a^. Our attention was first called 
to this inverted hierarchy by an observation that there are infrared divergences 
in a detailed study of the Kl vr^e+e" two-photon effect [33] which need to be 
canceled by the one- loop corrections to the radiative mode Kl — > 7r°7e+e~ through 
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the contributions of the soft radiative photons. This imphes that the theoretical and 
experimental analyses of —>■ 7r°e"'"e~ and Kl —>■ ■jTT^e'^e~ are tied together. The 
soft and coUinear photon regions of — > ^Tr^e'^e" form potential backgrounds to 
the studies of CP violation in the Kl 7r°e"'"e~ mode. 

The Kl —>■ 7r°7e"'"e~ mode also has an interest of its own. In recent years 
there have been important phenomenological studies of Kl — > 7r^77 in connection 
with ChPTh. This decay is calculable at one-loop (i.e., order E'^) ChPTh with no 
free parameters, yielding a very distinctive spectrum and a definite rate [34, 35]. 
Surprisingly, when the experiment was performed the spectrum was confirmed while 
the measured rate was more than a factor of 2 larger than predicted. The way 
out of this problem appears to have been provided by Cohen, Ecker and Pich 
(CEP) [36]. By adding an adjustable new effect at order E^, as well as including 
known corrections to the Kl tttttt vertex, they found that the predicted rate 
can be increased dramatically without modifying the shape of the spectrum much. 
This is also a surprising result, yet as far as we know it is the unique solution 
to the experimental puzzle. The ingredients of the mode studied in this thesis, 
Kl — > 7r°7e+e~, are the same as for Kl — > 7r'^77, except that one of the photons 
is off shell. Within the framework of the CEP calculation, the ingredients enter 
with different relative weights for off-shell photons. This will allow us to test the 
consistency of the theoretical resolution proposed for Kl — > 7r'^77. 

One of the goals of the next generation of rare kaon decay experiments is to 
attempt to observe CP violation in the decay Kl n°e'^e~. This reaction is special 
because we expect that direct CP violation (as opposed to the "mass matrix" CP 
violation already observed in the parameter e) may be the dominant component 
of the amplitude. This is in contrast with Kl — > tttt, where the direct effect is at 
most a few parts in a thousand of e. Direct CP violation distinguishes the Standard 
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Model from "superweak"-type models [37]. Moreover, the magnitude of the direct 
CP violation for this reaction is a precise prediction of the Standard Model, with 
very little hadronic uncertainty. In this thesis, we will update the analysis of the 
reaction Tr^e'^e", attempting to understand if we can be certain that an 

experimental measurement is in fact a signal of direct CP violation. 

One difficulty is that there are three possible components to the decay ampli- 
tude: 1) A CP-conscrving process which proceeds through two-photon exchanges, 
2) a mass matrix CP- violating effect proportional to the known parameter e, and 3) 
the direct AS — 1 CP-violating effect that we would hke to observe. The existence 
of the first of these shows that simply observing the total decay rate is not sufficient 
to unambiguously indicate the existence of CP violation. We need to either observe 
a truly CP-odd decay asymmetry, or else be confident on the basis of a theoretical 
calculation that the CP-conserving effect is safely smaller than the experimental 
signal. Unfortunately, the predictions in the literature for each of the components 
listed above exhibits a range of values, including some estimates where all three 
are similar in magnitude. However, the quality of the theoretical treatment can 
improve with time, effort and further experimental input. We will try to assess the 
present and future uncertainties in the theoretical analysis. 

There remain significant experimental difficulties before it is possible to mount a 
search sensitive to a branching ratio of a few times 10~^^. We will assume that such 
a sensitivity is reached. At the same time, it is reasonable to assume that we will 
have improved experimental information on the related rate Kl — > 7r°77, and that 
theoretical methods have provided a consistent phenomenology of this reaction. 
The CKM parameters will be somewhat more fully constrained in the future, but 
hadronic matrix element uncertainties will prevent a precise determination of the 
parameters relevant for CP violation, at least until B meson CP violation has been 
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extensively explored. With these expectations as our framework, will we be able 
to prove that the future experimental observation indicates the presence of direct 
CP violation? 

Our analysis shows that one will not be able to prove the existence of direct 
CP violation from the branching ratio for 7r^e~^e~ unless the decay rate for 

the related decay Ks 7r°e+e~ is also observed experimentally. This is yet more 
difficult than measuring the Kl decay, and poses a problem for the program of 
finding direct CP violation. It is possible but not certain, that the electron charge 
asymmetry can resolve this issue and, when combined with the rate, signal direct 
CP violation. 



2.2 Overview of the Analysis 

There is an extensive analysis associated with each of the three components of 
the decay amplitude for Kl vr^e+e" that were listed in the introduction, but 
only one work has appeared so far to treat the decay 7r°7e+e~ [38]. We will 

devote separate chapters of this thesis to each of these major issues. The purpose 
of the present chapter is to highlight the main issues that are to be discussed more 
fully later, and to indicate how they fit together in an overall description of the 
decay processes. 

For the first process, the direct CP component is, in a way, the simplest. The 
uncertainties are only in the basic parameters of the Standard Model, i.e., the mass 
of the top quark and the CKM parameters. The relevant hadronic matrix element 
is reliably known. Unfortunately the extraction of CKM elements has significant 

uncertainties, so that only a range of possible values can be given. This range 
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corresponds to Kl — > Tr'^e+e" branching ratios of a few times 10~^^. We discuss 
this range in chapter 3. 

The contribution of mass-matrix CP violation is more uncertain. The rate due 
to this source is given by the parameter e times the rate for Ks Tr^e'^e", so that 
the issue is the prediction of the CP-conserving Ks partial rate. Here the primary 
tool is ChPTh, with the pioneering treatment given by Ecker, Pich and de Rafael 
(EPR) [34, 35]. In chapter 4, we update their analysis, under essentially the same 
assumptions. The main new ingredient is the inclusion of the results of the one 
loop analysis of nonleptonic decays, which decreases the overall strength of the weak 
K ^ n transition. This yields a change in the weak counterterms and a decrease 
of the rate. However, more important is an assessment of the uncertainties of such 
a calculation, which we describe in chapter 5. Any such calculation has a range of 
uncertainties, most of which we are able to estimate based on past experience with 
chiral calculations. We systematically discuss these. Unfortunately we find that one 
issue in particular has a devastating sensitivity on this mode. In their analysis, EPR 
made an assumption that lies outside of ChPTh, that a certain weak counterterm 
satisfies W2 — 4:Lg where Lg is a known coefficient in the QCD effective chiral 
lagrangian. This assumption has no rigor, and the decay rate is very sensitive to 
the deviation W2 — 4:Lg. For any reasonable value of direct CP violation, there is an 
equally reasonable value of W2 that can reproduce the corresponding Kl — > 7r°e"^e~ 
decay rate. Given a measurement, we will then be intrinsically unable to decide if 
it is evidence of a nonzero value of direct CP violation or merely measures a value 
for W2- It is this what shows a need to measure the rate Ks — > 7r°e+e~. 

The third component is the CP-conserving amplitude that proceeds through 
the two-photon intermediate state Kl — > 7r°77, 77 — > e'^e", described in chapter 6. 
Here we must first understand the process Kl — > 7r°77. This has been calculated 
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in ChPTh at one loop order and has been measured experimentally. While the 
shape agrees with the chiral calculation, the theoretical rate misses by a factor of 
three. This has prompted some reanalyses of the theory of — > 7r°77, which we 
will take account of. However, the field has not reached a satisfactory conclusion 
on this mode, and it is clear that in the future the experimental and phenomeno- 
logical status of this reaction will undoubtedly improve. We study how possible 
resolutions of these analyses will influence the Kl 7r°77 decay rate. Ultimately 
this component should be satisfactorily understood. 

The ultimate problem is then our inability to distinguish, in a measurement 
of the Kl — > 7r*^e"'"e~ decay rate, the direct CP violation from the mass matrix 
effect. It is possible that the electron energy asymmetry may allow us to make 
this separation. The electron asymmetry comes from the interference of the CP- 
conserving two- photon process (even under the interchange of e+e") and the CP- 
violating one- photon process (odd under the e+e~ interchange). For many values 
of the presently favored parameter range, this asymmetry is very large i.e., of order 
50%. In this case its measurement is not far more difficult than a good measurement 
of the rate. If we in fact are able to reach an understanding of the two-photon 
process, through future phenomenology and experiments on n^'j^, then the 

asymmetry depends most critically on the CP-violating amplitude. If there is no 
direct CP violation, there is then a correlation between the Kl — > 7r°e+e~ decay 
rate and the electron asymmetry, parameterized by the unknown coefficient W2- 
As we detail in chapter 7, the presence of direct CP violation would upset this 
correlation, and in many cases would lead to a drastically different relative size 
of the asymmetry vs. decay rate, often even changing the sign of the asymmetry. 
Thus the asymmetry may be used to signal direct CP violation. Unfortunately this 
method is not foolproof. There exist combinations of values of W2 and CKM angles 
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for which the distinction between direct and mass matrix CP violation is not so 
great. In this case the analysis will be muddied by the inherent uncertainties in the 
theory. In chapter 7 we also explore the use of — Kg interference to sort out the 
direct CP violating amphtude. We outline the computation for the C(E'') ChPTh 
contributions to the process Kl —>■ n^je'^e' in chapter 8, and then we extend it to 
0{E^) in chapter 9. Finally, we recapitulate all our conclusions in chapter 10. 



CHAPTER 3 



DIRECT CP VIOLATION 



Direct AS = 1 CP violation is manifested in the "penguin" reactions pictured 
in Fig. 4. The QCD short distance corrections to this mode have been thoroughly 

analyzed to next-to-leading order by Buras et al. [39] (see also Ciuchini et al. 

[40]), and we will use their results. The primary weak operator responsible for the 
transition has the form 

T^W ~ = [C7v{lJ')Q7V + CtaQtaI , (3.1) 

where 



Qrv = {sd)v-A{ee)v, 

Q7A = {sd)v-A{ee)A- (3.2) 

The dominant contribution to the imaginary part of the coefficient Cn comes from 
the top quark, so that this is truly a short distance process. The coefficients have 
a CP-violating component 

Im Cr, = -Im (^^) y7i: (3.3) 

\ ud us / 
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Figure 4: "Penguin" diagram inducing AS = 1 CP violation. 



with the results of Ref. [39] yielding 



— = 0.743, — = -0.736, = 175, 

a a 



(3.4) 



with very httle dependence on Kjjg (the above is for k-jj^ = 0.3 GeV) and a 
negligible dependence on the low energy scale fi. 

The matrix element involved is well known via isospin symmetry from the 
charged current kaon decay, i.e.. 



(7r°(pJ|d7.«|X°(p^)) = ^(P^ + ^(Px -pJm, (3.5) 



V2 



V2 



with 



A = (0.65 ± 0.005) fm^. (3.6) 

The form factor /_ does not contribute significantly to the decay because its effect 
is proportional to m^. The decay rate is 
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BR{Kl ^ 7r°e+e-)dir = 4.16(ImAi)'(|/?A + Vrv), 
ImXt = ImVtaV* = | Kb 1 1 Kb I sin (5 = A^X\ 



(3.7) 



where K^ = I'^ubl sin (5, and A, A, r] refers to the Wolfenstein parametrization of the 
CKM matrix [9]. This results in 



BR{Kl 7r°e+e- 



jdir 



2.75 X 10-^2 X 



ImAi 
10-4 



for mt = 175 GeV. (3.8) 



The dependence on the top quark mass will of course be removed by a convincing 
precise measurement of mj. 

The CKM parameter Kb has the most favored values (including the recent 
CLEO data) [41] 



v: 



cb 



(3.9) 



0.0375 ± 0.0015 ± 0.0012 HQET [42] , 
0.036 ± 0.002 ± 0.003 B DHu models, 

0.039 ± 0.001 ± 0.004 B Xiu. 

The element Kb is nieasured by the inclusive decay B — > XuSiy in the electron 
endpoint region. The two inclusive calculations available yield 



ub 



K 



cb 



(3.10) 



0.082 ±0.006 [ACCMM] [43], 
0.074 ±0.007 [RDB] [44]. 
Models that calculate a set of exclusive decays {B — > Meu) can only be used to 
provide an upper bound on Kb since there are many final states (such as S — > ttttbi^ 
with TTTT nonresonant) that are not calculated. These limits are 



26 



ub 



0.12 [ISGW] [45], 
0.087 [BSW] [46], 
0.067 [KS] [47]. 
We will use the former measurements to estimate 



cb 



(3.11) 



ub 



cb 



= 0.078 ±0.007 ±0.010, (3.12) 

with the first uncertainty experimental and the second theoretical. In the Wolfen- 
stein parametrization of the CKM matrix, the values of V^^ and V^^ imply 



0.74 ±0.06, 



+ ^2 ^ 0.355 ±0.056. 



(3.13) 



Without any further analysis, these measurements imply an upper bound on ImA^: 



ImAt = |T/J]y^(,|sin(5 



< (1.0±0.3) X 10-Sin5. 



(3.14) 



A lower bound on this parameter can be found by consideration of the analysis of 
e. In the Wolfenstein parametrizations one has the approximate form 



|e| = (3.4 X 10-^)A^r]BK 



l + 1.3A\l-p) 



rut \ 



1.6 



mw 



(3.15) 



where Bk parametrizes the hadronic matrix element and is estimated to be in the 
range 0.33 < < 1. Using (1 — p) < 1.4 and B^ < 1 one finds for rrit — 175 GeV 



A^ri > 0.13, 



(3.16) 
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so that 

ImAt = A^X^ri > 6.8 x 10"^ (3.17) 

This brackets the range 

0.68 X 10"^ < ImAt < 1.3 X 10"^ (3.18) 

Note that ImAt is positive. These constraints yield a decay rate from direct CP 
violation of magnitude 

1.25 X 10-^^ < BR{Kl 7r°e+e-)dir < 4.6 x 10"^^ for = 175. (3.19) 
Alternately, the "best" values 

ImXt = 1.0 X 10-^ rrit = 175 GeV, (3.20) 
which we will take as our standard reference values, lead to a rate 

BR(Kl 7r°e+e-)dir = 2.32 x 10"-^^. (3.21) 

A more detailed analysis including a correlation between p and r) inherent in Eq. 
3.15, as well as the use of -B°S^ mixing (which constrains A^J {p — 1)^ + 77^ as 
a function of Jb) narrows the range only slightly because hadronic uncertainties 
dominate. 



CHAPTER 4 



REVISING THE EPR ANALYSIS 



In this chapter, we review the formahsm for analyzing mass matrix CP violation, 
first set forth by Ecker, Pich, de Rafael (EPR) [34, 35]. This amounts to the 
prediction of the decay rate for Kg — > 7r°e"'"e~, since the mass matrix effect is 
defined by 



^(XL^7r°e+e-)|MM = eA(i^5 ^ 7r°e+e-), 

e = (2.258 X 10-^)e^^/^. (4.1) 

We then redo the results taking into account recent work on the nonleptonic kaon 
decays to one-loop order. While this produces a significant numerical change, it 
is more important as a prelude to our subsequent analysis of uncertainties in the 
analysis. 

The prediction of Ks 7r°e"'"e~ comes from a comparison with vr+e+e", 
which contains many of the same ingredients. The reactions are displayed schemati- 
cally in Figs. 5,6. In these diagrams the round circles represent the electromagnetic 
coupling while the square boxes indicate the action of the weak interaction. We 
know the electromagnetic interactions of pions and kaons from direct measurement. 
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(c) 



Figure 5: Diagrams contributing to the reaction 7r+e+e^. 

Round circles represent the electromagnetic coupling while 
the squ£ire boxes indicate the action of the weak interac- 
tion. 

The weak X — > tt transition of Fig. 5a,b is known within some theoretical uncer- 
tainty from the use of chiral symmetry to relate it to X ^ 27r and K Stt. 
However, the weak Kttj vertex is not known a priori and needs to be extracted 
from the analysis of K'^ 7r+e+e~. 

The nonleptonic weak interactions are described by effective chiral lagrangians, 
organized in an expansion in powers of the energy, or equivalently in numbers of 
derivatives and masses. At lowest order, called order the physical transitions 
are described by a unique lagrangian 



C = GsTriX^D^UD'^U^), 
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Figure 6: Same as in Fig. 5 for Kl Tr^e+e . 

U = expfi— ^j, A^1...8, (4.2) 

where 0^ are the octet of pseudoscalar mesons (tt, K, rj). At next order, order 
the number of possible forms of lagrangians is quite large, and has been categorized 
by Kambor, Missimer and Wyler [48]. Not all of these contribute to X — > 27r 
and X — > TT, but certain linear combinations do influence these amplitudes. The 
formalism of ChPTh dictates that when an analysis is carried out to order E'^, that 
one use Eq. 4.2 at tree level, in which case one obtains from the X ^ tt decay rate 



^8 — ;^lKdK'j5'8) 
gT = 5.1. (4.3) 

Note that we neglect the CP violation in the nonleptonic amplitude, contained 
in Gg, because this is bounded to be tiny by the smallness of e'/e. 
In contrast, when evaluated at order E^, one must include one- loop diagrams in 
addition to the possible order E^ lagrangian. The loop diagrams involving tttt 
rescattering in the 7 = channel, K {1111)1=0 — ^ (7r7r)/=o as pictured in Fig. 
7, are quite large and are the major part of the order E^ analysis. While there 
is some ambiguity in the extraction of (see below), the enhancement from tttt 
rescattering leads to a smaller value of g^, with a good estimate being [49] 




Figure 7: One-loop diagrams involving tttt rescattering in the / = 
channel for Ks — > (7r7r)/=o — ^ (7r7r)/=o- The notation for 
vertices is as in Fig. 5. 



g'r = 4.3. (4.4) 

The X — > TT amplitude used in Figs. 5, 6 does not have the enhancement from tttt 
rescattering, and is given in terms of by 



^(K°^7r°) = -V2Gsk\ (4.5) 

We will explore further uncertainties in the K ^ ir vertex in the next chapter. 

To complete the diagrams of Figs. 5, 6 requires the electromagnetic vertices of 
kaons and pions. In ChPTh to order E'^ these are given by 



(7r+|J^|7r+) 
{K+\J''\K+) 



1 + r 



1 + r 



'2LM 
'2LM 




967r2F2 



(PK+PnT- 

(4.6) 



The first of these is known more fully from experiment, and has the form 
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(tt J^ tt ) = ^-r^' 

1 — /m^ 

m = 730 MeV. (4.7) 

Taylor expanding the latter form one determines Lq{^ — rrin) — (7.4 ± 0.7) x 10~^. 
The experimental charge radii 



(r2)^+ = (0.44 ± 0.02) fm^ 

(r2)^+ = (0.34 ± 0.05) W, (4.8) 

are compatible with this value. The final ingredient required for Figs. 5c, 6 is 
the weak photonic coupling. This includes both short distance and long distance 
physics, as illustrated in Fig. 8. While the short distance components have a 
reliable hadronic matrix element (it is due to the real parts of the coefficients 
discussed in the previous chapter), the QCD coefficient depends strongly on the 
low energy cutoff /i. In the full matrix element, this dependence is canceled by a 
corresponding dependence on /i of the long distance physics. However, since the 
long distance physics is not calculable, we must attempt to determine this coupling 
phenomenologically. The innovation of EPR was to elucidate the possible forms 
that this coupling could take. They found that there were two possible chiral 
lagrangians that could contribute to this process: 

Cw = '-^F^''[w^Tr{QXe-^7C^C,)+W2Tr{Q/:^Xe-i7C,)], 

= -{d^U-te[A^,U])U^ = -{D^U)Ul (4.9) 

In the presence of the short distance electroweak penguin effect due to Z exchange 
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Figure 8: Diagrams contributing to the short distance weak photonic 
coupling for K — > 7re+e". 

with an axial electron coupling, we need a third effective lagrangian, not present in 
EPR, 



C^w = -^G^w^e^^^^eTr . (4.10) 

The correspondence with the notation of the previous chapter is 



Im W5 = ^ = ^ImA*. (4.11) 

47r|T/,T/*|^3 a 

Note that the labeling of takes into account the lagrangians labeled by ^3, ^4 
defined by EPR [34, 35, 50] that contributed to other radiative K decays. In fact 
this form of £J^ is closely related to w\ term since by using identities of the U 
matrix it can be shown that 



= {-e)e-fi^eTr{QXe_i7C^) 
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= Q)e7^eTr(A6-i7>C.), (4.12) 



where in the second hne we have integrated by parts thus subsequently used the 
equation of motion so that 

^-^^F^^WiTriQXe^ijCi^C^) = wie^u,eTr{\Q-„Cy). (4.13) 

This allows us to identify the short distance CP violating part of wi. 

Im = ^ ^ ^ImAt. (4.14) 

The real parts of Wi, W2 contain long distance contributions and hence are not 
predictable by present techniques. They need to be extracted from experimental 
measurements. The other process available for this procedure is 7r+e+e~ 
(unless Ks —>■ Tr^e+e" is measured in the future). Unfortunately one cannot fix 
both wi, W2 in this way, so that one must add a theoretical assumption in order to 
proceed. The relevant amplitudes are 
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M{Ks ^ TT^'e'^e ) = — —ds{pK + Pj^uji,v, 

MiKL^TT^e+e-) = -^{Pk + P.T'aldvj^ + dAj^j5]v, (4.15) 

with 



d+ = «;+ + 0ftr(5^) +0^(g^), 
ds = Rews + 2(f)K{q^), 

w+ = — [wl + - 12LI) - g log ^4 , 



35 



1 TD 

ws = w+ + 167r2K-4L^) + -log^, 



''K 



dx 



— 7jX(l — x) 



log 



— T^xil — x) 
mi 



(4.16) 



and for the CP violating decay 



dv — eds ilm wi, 

3 

dA = — — ilmw5. (4.17) 

The goal of the search for direct CP violation is to separate the Im wi,5 terms from 
the mass matrix effect eds- Note that in these expressions we have neglected the 
possible direct CP violation in the X — > tt transition (which is bounded to be very 
small by the measurement of e'/ e) and the contribution of Re to CP-conserving 
decays (since Re ^5 ~ Im ^5 << Re wi). 

The EPR analysis of K'^ 7r+e+e~ uses the tree level value of (^g, gl^^^ = 5.1. 
The decay rate is consistent with two values of Re and a subsequent analysis 
of the decay spectrum favored the lower value for Re twi, i.e. Re ty+ = 1.16 ± 0.08 
[34]. However, given that one is working to one-loop order, it is more consistent to 
use the one-loop value for g^, g^^""^ = 4.3. Because of the presence of the Lg term, 
this is not just a rescaling of the value of w+. An additional change that we make 
is to use the known full electromagnetic vertex in the pole diagrams, rather than 
just the first term in the expansion of the form factor. Note that because of the 
factor of -p^ in the weak matrix element, the only significant form factor is that 
of the pion in Fig. 5. This implies the replacement 



^ F2 (1 - q^ml) 



(4.18) 
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in the formula for The associated logarithm with Ri are much smaller than 
the Lg dependence and are not influenced much by this replacement. As a technical 
note, we comment that some potential modifications using a phenomenological pion 
form factor could lead to a lack of gauge invariance. By modifying the coefficient 
of a gauge invariant effective lagrangian, we preserve the gauge invariant nature of 
the amphtude. With these changes, we find 



[Without the second change, we would have had Re = 1.33 ± 0.065]. This is 
illustrated in Fig. 9. 

One cannot simply transfer this information to Ks or Kl decays, because a 
different linear combination enters 



However, EPR deal with this problem by making the assumption that W2 — 4:Lg 
resulting in Re = 0.73 ± 0.08 [34]. They note that this assumption is not part 
of ChPTh, but do not explore the consequences if it is not correct. We will discuss 
this in the next chapter, finding a very strong sensitivity. At this stage we note 
that if one makes the assumption of W2 — 4L9, one obtains Re Wg = 0.58 ± 0.10 for 
our value of Re w+. 

At this value oi Wg, the direct and mass matrix contributions are comparable: 



Re w+ = 1.01 ±0.10. 



(4.19) 




(4.20) 



dv — eds — 



ilm Wi 



3 




dA = il.O X 10 



,-3 



(4.21) 
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Figure 9: The branching ratio BR{K'^ 7r+e+e~) is plotted against 
Re The soUd curves £ire obtained using the extreme 
values of the error intervals of \Vud\ = 0.9753 ± 0.0006, \Vus\ 
= 0.221 ± 0.003 and Lg{m,i) = (7.4 ± 0.7) x 10-^ while the 
dashed curve corresponds to the central values. The ex- 
perimental value of the branching ratio ± its experimental 
errors are indicated by deished horizontal lines. 

when evaluated with rrit — 175 GeV, ImA^ = 10~^. This leads to a branching ratio 



BR{Kl 7r°e+e-)MM = 0.37 x ID" 



-12 



(4.22) 



if there is no direct CP violation (ImAj = 0) [EPR found BRmm = 1-5 x 10 in 
this case [35]], vs. 



BR{Kl 7r°e+e-)cp = 1.78 x 10" 



12 



(4.23) 
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Figure 10: The branching ratio BR{Kl — > 7r°e+e~)MM is plotted 
against Re The solid curves are obtained using the 
extreme values of the error intervals of \Vud\ = 0.9753 ± 
0.0006, iKsl = 0.221 ± 0.003 and Lg{m,i) = (7.4 ± 0.7) x 
10~^, while the dashed curve corresponds to the central 
values. The experimental limit on Re w+ ± its experimen- 
tal errors are indicated by deished vertical lines. 

for the full set of parameters given above in Eq. 4.16. Figs. 10 and 11 show the 
dependence of both the above branching ratios on Re w+. The addition of mass 
matrix CP violation in this analysis led to a small decrease in the rate compared 
to the purely direct CP violation of the previous chapter, Eq. 3.21, because of the 
cancellation in the imaginary part of dy- However, this may change if W2 7^ 4L9. 

One of the results that we will see in the next chapter is that the mass matrix 
contribution to the branching ratio is near a minimum value when W2 is close to 



39 



BR(K2^7T°e^e")cpX10 




0.0 0.5 1.0 1.5 

RevrJ 



Figure 11: Same as Fig. 10 for the branching ratio BR{K^ 7r°e+e )cp. 

4L9. For other values of the rate can easily be an order of magnitude larger. 
Although our value for the branching ratio is a factor of four below that of EPR, 
both estimates are similar in saying that the mass matrix contribution will be small 
as long as W2 = 4L9. 

Using the tree level value for g^, gf"^^ = 5.1, the mass matrix contribution to 
the branching ratio is 



BR{Kl 7r°e+e-)MM = 0.55 x lO'^^. (4.24) 



CHAPTER 5 



UNCERTAINTIES IN MASS MATRIX CP VIOLATION 

Our goal in this chapter is to assess how well we understand the prediction for 
Ks TT^e+e". The most important effect will be discussed in section 5.3 below, but 
we proceed systematically to discuss even contributions that have less uncertainty. 

5.1 Purely Electromagnetic Vertices 

The electromagnetic vertices enter in diagrams 5a,b, 6. The uncertainty here is 
in the choice of whether to use the chiral expansion of the form factor truncated 
at order q^, Eq. 4.6, or the full dependence of the monopole form factor, Eqs. 
4.7, 4.18. The first choice is natural when one is working to a given order in the 
chiral energy expansion, but the latter choice clearly includes more of the physics 
that is known about the electromagnetic vertex. Note that essentially only the 
pion form factor is relevant, because the diagram involving the kaon form factor. 
Fig. 5b, is suppressed by a factor of m^/ml- with respect to Fig. 5a because of the 
momentum dependence of the weak K ^ tt transition. The use of the full form 
factor produces a modest variation in the value of w+ (i.e. Re w+ — 1.01 instead 
of Re w+ — 1.33). Because of cancellations in the amphtude this provokes a 
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more extreme variation on the decay rate. The lowest order chiral vertex, Eq. 4.6, 
produces a decay rate 

BR{Kl 7r°e+e")MM = 1-96 x 10"^^ (5.1) 

instead of the result of Eq. 4.22. [We note that if we had also modified W2 in 
the same way as 4Lg as in Eq. 4.18 we would have a BR{Ki — > 7r°e"'"e~)MM — 
1.33 X 10~^^.] While we feel that it is good physics to use the full electromagnetic 
form factor, one could also interpret these results as an uncertainty in the analysis 
due to higher order terms in g^, with that uncertainty being of order 2 x lO"-*^^. 



5.2 The Weak K ^ Vertex 

We have already given one indication of the sensitivity of the result to the size of 
the K ^ TT transition. Under otherwise identical assumptions, g^'^'^ = 5.1 yielded 
the rate in Eq. 4.24, while g^^°^ — 4.3 produced the result of Eq. 4.22. These 
modifications to also include corresponding modifications to the Ktt^ vertex 
required by gauge invariance. This is automatically maintained, however, by the 
use of gauge invariant effective lagrangians. In order to appreciate that this change 
in g^ is not the only uncertainties in the Ktt amplitude, one needs to understand a 
bit more about the chiral phenomenology of — > 27r and K ^ Stt. 

Chiral symmetry relates processes with different numbers of pions, such as X — > 
TT vs. X — > Stt. The predictions are compactly contained in the chiral lagrangians, 
but can also be obtained using the soft pion theorems, which was the methodology 
used in the 1960's. The only advantage of the latter technique is that it relies only 
on chiral SU(2) while modern chiral lagrangian analyses have always involved chiral 
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SU(3) symmetry. [Presumably the latter could be reformulated in chiral SU(2), but 
no one has yet done this.] The soft pion analysis indicates that one obtains the 
same relation (up to terms of order m^) between X — > 27r and X — > Stt for any 
lagrangian that survives in any soft pion limit oi K ^ Sir (i.e., pi 0). The 
only lagrangians that do not survive in any soft pion limit involve four separate 
derivatives on the four fields of K — > Stt, e.g. 

Auaxtic = j^Tr (XeD^UDM^D'^UD'^U^) . (5.2) 

This lagrangian yields a matrix element proportional to {px ■ Pi){p2 ■ Ps) which 
clearly vanishes as any pj — > 0. In contrast most of the order lagrangians do 
not vanish in all soft pion limits. An example is 

£' = j^Tr (U^ XeDi^D^D^U^ D'^U) , (5.3) 

which then yields the same relations of X ^ Stt and ^ 27r as does the lowest 
order result C given in Eq. 4.2. [This phenomenon is explained in more detail in 
Rcf. [51].] Since the only inputs to the chiral phenomenology arc the amplitude for 
K ^ 27: and K Stt, it follows that one cannot distinguish a combination C-\- C 
from a lagrangian involving C only. 

However, when we discuss the K ^ n vertex, there is a distinction between 
these various lagrangians. For example C in Eq. 5.3 involves a minimum of three 
meson fields, and hence contributes to X ^ 27r and X — > Stt but not at all to 
K ^ TT. In contrast, the lowest order lagrangian Eq. 4.2 contributes to all of — > 
TT, X — > 27r, K — > Stt. Since known phenomenology cannot distinguish between 
linear combinations oi C + C, this manifests itself in an uncertainty in the K ^ tt 
vertex. Since the higher order lagrangians are known to make a 25% difference 
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in relation between X — > 27r and K — > Sir, it would be unreasonable to take this 
uncertainty in — > tt to be any less than 25-30%. Using gg = 4.3 x (1 ± 30%) 
yields a range 

BR{Kl 7r°e+e-)MM ^ 0(10"^^) ^ 1.5 x lO'^l (5.4) 

This again indicates that the analysis has significant cancellations present and that 
modest variations in the analysis can lead to uncertainties of order 1.5x10"^^. 

5.3 The Kttj Vertex 

The parameter w+ was determined from the analysis of K'^ — > Tr'^e'^e" . This 
arrangement has some intrinsic uncertainty because it was performed to a given 
order in the chiral energy expansion. It would be reasonable to take this uncertainty 
at 30%. 

However, this uncertainty is dwarfed by the errors introduced by the assumption 
of W2 — 4L9. There is nothing about chiral symmetry that forces such a relation. 
For example, if the process of Fig. 12 were to contribute to the weak coupling, 
the relation W2 = 4L9 would not occur except for a special value of the K ^ ai 
amphtude. This is highly unlikely, and we can easily accept W2 — 4L9 ^ 0. The 
crucial distinction here is between models and rigorous theory. ChPTh is a rigorous 
method which expresses true relationships in QCD to a given order in the energy 
expansion. However, an assumption such as ^2 = 4L9 may be true or false in a way 
that we cannot decide based on QCD. It may occur within some models, yet we 
have no guidance as to whether it is correct in nature. We cannot base something 
as important as the observation of direct CP violation on something as fiimsy as a 
model. 
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Figure 12: Diagram contributing to the wecik photon coupling. The 
notation for the vertices is as in Fig. 5. 

Unfortunately the decay rate for Ks — > Tr^e'^e" depends very strongly on the 
value of W2- This strong dependence was observed by Littenberg and Valencia 
[21]. If we were to chose W2 — 0, the rate would be two orders of magnitude 
larger. In Fig. 13 we plot the branching ratio for the mass matrix contribution to 
Kl — >■ TT^e'^e" vs. W2. We see that reasonable values of W2, we get a wide range of 
values of the branching ratio. Conversely, a measured value of BR{Kl Tr^e+e") 
in the range of 10~^^ few x 10~^^ could be interpreted in terms of a reasonable 
value of W2- This is then an enormous uncertainty in the mass matrix contributions 
to Kl n^e+e-. 

This uncertainty could be removed if one measured the rate of Ks vr^e+e". 
The mass matrix contribution would then be known, see Eq. 4.1. However, this 
task is not easy experimentally. 
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Figure 13: The mass matrix contribution for the branching ratio 
BR{Kl — > 7r°e"'"e~)MM is plotted against W2- The convention 
for the soUd and dashed curves is the same as for Fig. 9. 



CHAPTER 6 



THE CP-CONSERVING AMPLITUDE 

The Kl 7r°e+e~ transition can also take place through a CP-conserving two- 
photon intermediate state. If we ignore the electron mass, the form of the amplitude 
will be 

M{Kl ^ 7r0e+e-)cpc = G^a^Kp^ ■ {k - k'){pj, + p^rn^^v, (6.1) 

where X is a form factor and the extra antisymmetry under k ^ k' {k = ke+, k' 
= kg-) is a reflection of the properties under a CP transformation. In order to 
calculate this we need to understand the Kl — > 7r°77 transition first. 

We are fortunate that — >■ 7r°77 is accessible to present experiments, and 
significant phenomenology has been performed on this reaction. We will utilize 
the work of Cohen, Ecker and Pich [36] as representative of present work, with 
the understanding that future experimental and theoretical work will clarify the 
analysis considerably. 

The most general form of the K 7r77 amplitude depends on four independent 
invariant amplitudes [35] A, B, C and D: 
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+2— (pk ■ Qi (I2PK +Pk-(12 QiPk - Qi ' IiPkPk -Pr-QiPk' ?2 9'"') 

+Ce'^'"'''q,^q,, (6.2) 
^u,upa . +pj^-q^ q^p) Pk^ 



D 



In the limit where CP is conserved, the amphtudes A and S contribute to 
K2 — > 7r°77 whereas Ki — > 77^77 involves the other two amplitudes C and D. All 
four amplitudes contribute to — > 7r"'"77. Only A and C are non-vanishing to 
lowest non-trivial order, in ChPTh. 

When the photons couple to e+e~, as in Fig. 14, it is well known that the 
A amplitude contributes to 7r°e"'"e~ only proportionally to mg, which is a 

small effect that we will drop. It is the B amplitude that is important for the e^e" 
final state. The authors use a representation that fits the known Kl — > 7r+7r~7r'^ 
amplitude in a dispersive treatment of — > 77^77 and find 
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Figure 14: CP-nonviolating diagram involving two photons coupling 



The most important ingredient above is ay which is a parameter representing 
the vector meson exchange contributions to the A, B and D amplitudes. A fit 
(shown in Fig. 15) to the decay rate [using BR[Kl — > 7r°77) = 1.7 x 10~^ [52]] 
and the 77 spectrum in Kl 71^77 (depicted in Fig. 16) indicates a value around 
tty — —0.96. This parameter was very important in increasing the chiral prediction 
of the decay rate to be in agreement with experiment. We have explored the 
sensitivity of this parameter to changes in the analysis and have found that 25% 
changes in the dispersive treatment lead to a factor of 2 change in a^, so that this 
value is still quite uncertain. 

Several authors [50, 36, 53, 54, 55] have calculated the contribution of the on- 
shell two-photon intermediate state to Kl 7i^e^e~ . Although this is sometimes 
referred to as the absorption contribution, it is not the full absorption part since 
there is a further cut due to on-shell pions. Besides this, the full CP-conserving 
amplitude also receives contribution form the dispersive part of the amplitude. 



to e+e . The notation for the vertices is as in Fig. 5. 




x> 1. 



(6.3) 
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Figure 15: 77-inv£Lriant distribution for Kl 7r°77: 0{E'^) (dotted 
curve), 0{E^) with ay = (dashed curve), 0{E^) with ay 
= —0.96 (full curve). The spectrum is normalized to the 
50 unambiguously events of NA31 [52] (without accep- 
tance corrections). 



with off-shell photons (and pions). The calculation of this is complicated by the 
sensitivity of the loop integral to high momentum, as the Feynman diagram of Fig. 
14 will diverge if we treat the B amplitude as a constant in qf and g|. However, 
the remedy to this is well known; the couplings of virtual photons to hadrons is 
governed by form factors which lead to suppression of the couplings at high q^. We 
will include an estimate of these form factors and this will allow us to calculate the 
dispersive component of the CP-conserving amplitude. 
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Figure 16: 77-invciriant mass distributions of the 50 unambiguously 
Kl — > 7r°77 events reconstructed by NA31 [52] (his- 
tograms). The dashed Hne shows the estimated back- 
ground. The dotted hne simulates the 0{E^) ChPTh 
prediction. Crosses indicate the experimental acceptance 
(scale on the right). 

The two-photon loop integral in the limit mg — > is given by 
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7r°e+e )cpc 
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(6.4) 
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Here F(g^) are the form factors for the virtual photon couphngs. The structure 
above is certainly an approximation, as in general the virtual photon dependence 
need not be only an overall factor of F{q^). However, the above form would be 
sufficient to capture the kinematic variation if only one photon is off-shell [given an 
appropriate F{q^)]. Since we only need a minor form factor suppression to tame 
the logarithmic divergences, we feel that this structure will be sufficient for our 
estimate. We choose 

Hq') = (6.5) 
- my 

which is a good representation of almost any mesonic form factor, with my ~ m^. 
Neglecting terms that are suppressed by powers of 1/m^, we find the amplitude of 
Eq. 6.1 above, with 



^ B{x) 



e , 



(6.6) 



167r2m|^ 

where s — {k-\-k')'^. The log factor is of course expected, since the photon "ab- 
sorptive" part comes from the expansion log(— ,s) = logs + m. 

This representation of the amplitude leads to a CP-conserving branching ratio 

of 



BR{Kl 7r°e+e-)cpc = 4.89 x 10"^^ (6.7) 

for ay = —0.96. More generally we show the CP-conserving branching ratio vs. 
ay in Fig. 17. Note that while for many values of ay the CP-conserving rate is 
small compared to the CP-violating rate of previous chapters, these two rates are 
comparable for some range of parameters. Note that there is no interference in 
the rate between the CP-conserving and violating components, so that the rates 
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Figure 17: The CP-conserving branching ratio BR{Kl — > 7r°e"'"e~)cpc 
is plotted against ay. The convention for the solid and 
dashed curves is the same as for Fig. 9. The assumed 
value for ay is indicated by the dashed vertical line. 



just add, as shown in Fig. 18. However, there can be a CP-odd asymmetry in the 
electron positron energies, to which we turn in the next chapter. 
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Figure 18: The complete branching ratio BR[Kl vr^e+e") is plot- 
ted against W2- The convention for the solid and dashed 
curves is the same as for Fig. 9. ay — —0.96 is assumed. 



CHAPTER 



7 



THE ELECTRON ENERGY ASYMMETRY AND 



TIME-DEPENDENT INTERFERENCE 



7.1 The Electron Energy Asymmetry 

When both CP- violating and CP-conserving amphtudes contribute to the decay, 
there will be an asymmetry in the electron-positron energy distribution 



This will be quite large if the two amplitudes are comparable. In contrast to the 
overall decay rate, this asymmetry is an unambiguous signal of CP violation. It 
may be even more useful if it can be used to prove the existence of direct CP 
violation. This can occur because the asymmetry is sensitive to the phase of the 
CP- violating amplitude, and mass matrix CP violation has a unique phase (that 
of e), while direct CP violation will in general have a different phase. 

For the electron energy asymmetry to be useful as a diagnostic of the form of 
CP violation, the CP-conserving two-photon amplitudes must be well known. As 
we discussed in the previous chapter, it is reasonable to expect that this will be 
true in the future after further phenomenology of 77^77. For illustrative 



N{E+ > E_) - N{E+ < E_) 
N{E+ > EJ) + N{E+ < E_) ' 



(7.1) 
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purposes, we will use Uy — —0.96. The analysis below will need to be redone in the 
future if this value of ay changes significantly. However, the pattern of the analysis 
and the general conclusions will be valid for a wide range of values of ay. 

The amplitudes involved have been given in the previous chapter. Note that 
the two-photon amplitude has both real and imaginary parts, the mass matrix 
CP- violating amplitude has a phase of 45°, while the direct component is purely 
imaginary. The asymmetry is proportional to the imaginary part of B times the 
real part of dy minus the product of the real part of B and imaginary part of dy. 
The asymmetry may be defined in differential form: 



— (E+>E) = /^^'^''''''^dy-^ 

dz Jo dydz ' 

t^r,^ ^ , /-o , dT 

— (E+ <E_) = , dy——, 
dz J-^\^{i,z,r2) dydz 

A(l,z,r^) = l + z'^ + r^-2z-2r^ -2r^z, 



A{z) 



rriK 

dT{E+ > E_)/dz - dT{E+ < E_)/dz 



dT{E+> E_)/dz + dr{E+< E_)/dz' 
z = {k + k'f/mj,. (7.2) 

In Fig. 19 we plot the differential asymmetry for several values of W2 in the 
case when there is no direct CP violation. Fig. 20 gives the same information for 
ImXt = 10~^. We see that the asymmetry is sizable for many values of W2 and that 
it depends significantly on direct CP violation, even possibly changing sign. The 
integrated asymmetry is plotted versus W2 in Figs. 21, 22. 

Both the decay rate and the asymmetry depend on Wg and this forms the main 
uncertainty in the analysis. However, it is possible to remove this uncertainty by 
measuring both observables. For a given value of the CP-conserving amplitude. 
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Figure 19: The differential asymmetry A{Kl n^e'^e') is plotted 
against z = {k + k'Y/m\ for ^2 = 1 xlO~^ (curve 1), 2 
xlO~^ (curve 2), . . ., 5 xlO~^ (curve 5), in the ceise when 
there is no direct CP violation, ay = —0.96 is assumed. 



there is a strict correlation between the two. In Fig. 23 we plot the values of the 
branching ratio and the integrated asymmetry A as one varies W2. We see that 
the curves with and without direct CP violation are well separated for much of the 
range. To the extent that we understand the CP-conserving amplitude, this can 
be used as a diagnostic test for direct CP violation. 
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Figure 20: Same as Fig. 19 for ImAj = 10 ^ 



7.2 Time-Dependent Interference of i^i, K2 — ^ tt e'''e 



Littenberg [56] has suggested a time-dependent analysis of a state starting out 
as a in order to extract maximal information about the decay mechanism. This 
is far more demanding experimentally then simply measuring Kl or Kg decays 
separately. However, we analyze this technique in order to assess its usefulness. 

A state that at i = is a K*^ will evolve into a mixture of and Kg'- 



\K%t)) = ^Y^^ {e-^^-* m) + e\K,)] + e^^^* [\K,) + e\K,)]} , (7.3) 
where 
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Figure 21: The integrated asymmetry is plotted vs. W2 in the case 
when there is no direct CP violation. The convention for 
the solid and deished curves is the same as for Fig. 9. 
ay — —0.96 is assumed. 



(7.4) 



Ignoring small effects such as second order CP violation and direct CP- violating 
components in the parameter e, we then have a time development proportional to 



1 

— 

2 ^ 

+ 2Re 



Ac 



(tAs + Adir + Acpc) A*se-'^^'^-^s)t 



e 2 



(7.5) 
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Figure 22: Same as Fig. 21 for ImAj = 10~^. 



Here A^u, eAs and ^cpc are the amplitudes analyzed in chapters 3, 4 and 6 
respectively. Measurements at early time {t <^ l/r^) determine T{Ks vr^e+e") 
while at the late time {t » l/r^) one observes r{KL — > 7r°e+e~). These contain the 
information described in preceding chapters. However, in the interference region t 
— O ^l/(m^ — irig)^ ~ 0{1/Tg), we obtain extra information about the separate 
contributions to the decay amplitude. 

In Fig. 24 we show the time-dependent signal for the cases of pure mass matrix 
CP violation and the addition of direct CP violation, using the analysis of the 
previous chapters with W2 — 4:Lg. We see that the shape of the interference region 
does differentiate these two cases, but that for the case studied the dependence on 
direct CP violation is not so large as to allow an easy experimental determination. 




Figure 23: The integrated asymmetry A{Kl — > Ti^e^e~) is plotted vs. 

the branching ratio BR{Kl vr^e+e") as one varies W2. 
ay — —0.96 is assumed. The soHd curve describes the ceise 
when there is no direct CP violation, the dashed curve is 
for ImAt = 10-^. 
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Figure 24: The normalized time distribution of — > 7r°e"'"e~ decays 
is plotted vs. time t in nsecs. oy = —0.96 is cissumed. The 
solid curve describes the case when there is no direct CP 
violation, the dashed curve is for ImAf = 10~^. 



CHAPTER 8 



THE 0{E^) CALCULATION 



First let us provide the straightforward O(E^) calculation oiA4{KL — > 7r^^e'^e~) 
within ChPTh. This is the generalization to kf ^ of the original chiral calculation 
of EPR [34, 35] . Here ki is the momentum of the off-shell photon. This captures 
all the k\/m1 and k\/m\ variations of the amplitudes at this order in the energy 
expansion. There can be further k\/{l GeV)^ corrections which correspond to 
O(E^) and higher. The easiest technique for this calculation uses the basis where 
the kaon and pion fields are transformed so that the propagators have no off- 
diagonal terms, as described in Ref. [34, 35]. The relevant diagrams are then 
shown in Fig. 25. Defining 'g as 

g = G8/3, G, = G^lK^dCl^s, l^sl ^ 5.1, (8.1) 

the diagrams give the following integrals, respectively: 

_.2^, ( dH ?,[{p^-p,f -ml]-2[{f -ml) + {l-k,-k,f -ml] 
gg,, J ^2^^^ _ _ _ _ 

(8.2) 



2- f dH 3[{p^-p^f-ml]-2[{l + k,y-ml+il-k2r-ml] 



(27r)4 (/2 - ml)[{l + kif - ml][{l - k^f - ml] 
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K,TT 




c) 



d) 



Figure 25: Diagrams relevant to the process — > 7r°7e"'"e at 0{E^) 
and 0{E^). 



.3) 



dH 1 



(8.4) 



+ 



(27r)4 [ (/2 - m2)[(/ _ fc^)2 _ ^2] _ ^2)[(; _ ^2)2 _ ^2 



(8.5) 



Interestingly when we add these together the K ^ Sti amplitude factors out from 
the remaining loop integral resulting in 
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dH [g,,{P - ml) - {21 + h),{2l - k^),] 



(27r4) (/2 - ml)[{l + A;i)2 - ml][{l - k2f - m: 



(8.6) 



It is not hard to verify that this result satisfies the constraints of gauge invariance 
kiM^v = k2M.^u = 0. At this stage, the integral may be parametrized and inte- 
grated using standard Feynman-diagram techniques. Let us keep photon number 
one as the off-shell photon and set A;| = 0. In this case the amplitude with one 
photon off-shell is described by 



Ml^ = Qe^-9[{PK-Pof - ml\ (^) ^l^^^^^l^^^M^^ + 2/«)], (8.7) 



with 





{F{s) - F{kl)\ 




{G{s)-G{kl)\. 




The notation is defined by 



s = {Px-Pof = {kl + k2f 



(8.9) 



and 




(8.10) 




1 



(8.11) 



The above functions are related to those presented by CEP [36]: 
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F{a) 



2ml 



CEP 



(8.12) 



G{a) 



2ml 



RCEP 



a \ 1 
4ml) ^ 6. 



remembering Eq. 6.3: 



1 - 


1 r 




X - 


1 + 


1 

Ax 


1 

~6 


+ 7 


1 

~6 


+ 7 



arcsm 



X<1, 



log , + tn 

1 + Jl-l/x 



1 — \Jl/x — 1 arcsin (v^) 



X > 1, 
X < 1, 



; , 1-Jl-l/a; 

1 + ^1 - l/x log ^, + m 



1 + Jl - 1/a; 



(8.13) 



x>l. 
(8.14) 



This agrees with the EPR result in the A;^ — > limit. 

At this order we have also calculated the additional contribution resulting from 
the kaons circulating in the loops of Fig. 25. They give rise to 

(8.15) 

The resulting integral is similar to that of Eq. 8.8, substituting the mass of the 
pion with that of the kaon. Attaching an e+e" couple to either photon and adding 
all the above contributions together, the result we obtain for the branching ratio is 



BR{Kl 7r%e+e-) = 1.0 x 10" 



(8.16) 



66 



dBR(KL^7TV^e")/dz 
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Figure 26: Differential branching ratio dBR{KL 7r°7e+e-)/(i2; to C>(E^). 



With the definitions 



— ' y 



Pk ■ (^1 - ^2) 



(8.17) 



the decay distributions in z and y provide more detailed information. We present 
them in Figs. 26 and 27. 
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CHAPTER 9 



THE 0{E^) CALCULATION 



We also wish to extend this calculation along the lines proposed by CEP [36] , 
who provide a plausible solution to the problem raised by the experimental rate 
not agreeing with the O(E^) calculation when both photons are on-shell. The two 
primary new ingredients involve known physics that surfaces at the next order in the 
energy expansion. The first involves the known quadratic energy variation of the 
X — > Stt amplitude, which occurs from higher order terms in the weak nonleptonic 
lagrangian [51, 55, 57] . While the full one-loop structure of this is known [48, 49, 58] , 
it involves complicated nonanalytic functions and we approximate the result at 
O(E^) by an analytic polynomial which provides a good description of the data 
throughout the physical region: 

M{K ^ 7r+7r"7r°) = 4aip^ • PoP+ • P- + 4a2(PK • P+Po ■P-+Pk- P-Po ' (9-1) 
using 

oi = 3.1 X and 02 = -1.26 x IQ-^m^l (9.2) 

ai and 02 are obtained from a fit to the amplitude for Kl — > 7r°7r"^7r~ [51] and 
to the amplitude and spectrum for — > 7r°e+e~ [36], so that their values are 
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constrained within their theoretical uncertainty of 10 - 20%. We have numerically 
verified that such a variation of said parameters involves a very modest change in 
the shape of the spectrum for — > 77r°e''"e~ and a change in its final branching 
ratio somewhat smaller than the uncertainty on the parameters. 

The other ingredient involves vector meson exchange such as in Fig. 28. Some 
of such contributions are known, but there are others such as those depicted in Fig. 
29 which have the same structure but an unknown strength, leaving the total result 
unknown. In Ref. [36] the result is parametrized by a "subtraction constant" that 
must be fit to the data. 

In principle one can add the ingredients to the amplitudes and perform a dis- 
persive calculation of the total transition matrix element. In practice it is simple to 
convert the problem into an effective field theory and and do a Feynman-diagram 
calculation which will yield the same result. We follow this latter course. 

The Feynman diagrams are the same as shown in Fig. 25, although the vertices 
are modified by the presence of 0(E'^) terms in the energy expansion. Not only 
does the direct K —> Sir vertex change to the form given in Eq. 9.1, but also the 
weak vertices with one and two photons have a related change. The easiest way to 
determine these is to write a gauge invariant effective lagrangian with coefficients 
adjusted to reproduce Eq. 9.1. We find 

Mi,(K 7r+7r"7r%) = 4aie(p+-p_)^p^-po+4a2e(p+-p_)^(PoPK^+PKPoM)' (9-3) 
M^^{K 7r+7r"7r%7) = Saie^g^^pj^ ■ + 8a2e^(p^^Po^ + PKuPof^)- (9-4) 




Figure 28: Vector meson exchange diagrams contributing to Kl — > 

The resulting calculation follows the same steps as described above, but is more 
involved and is not easy to present in a simple form. We have checked that our 
result is gauge invariant and reduces to that of CEP in the limit of on-shell photons. 

The contribution proportional to ai can be computed analogously to those al- 
ready calculated for the C(E^) case: 

= 4aie2(z - 2r^)(l + - z)^{g^,h ■ - k2^.h,)[l + 27(m^)], (9.5) 



a) 



V = p,u) 




b) 



Figure 29: Vector meson exchange diagrams contributing to Kl 
7r°7e+e~ with unknown strength. 



where 



mi<- mj^ mj( 



The a2 part originates another set of integrals which can be written as 



72 



dH { {2l + ki)^{2l-k2)u{l + h)p{l-k2)a 



+ 



{27rY [ (/2 - ml)[{l + ki^ - ml][{l - k2Y - ml] 
{21 + k2),{2l - ki)^{l + k2)p{l - ki)^ \ 



{P - ml)[{l - A;i)2 - ml][{l + k2f - ml] J 

(9.8) 



r dH 1 

M% = Sa2ipK^J)ou+PKuPof^y J (^27rYP-ml ' ^^'^^ 



, .rl . , a a ^ 1 i (21- kl) J2l - kl)a 

- ^a2{PoPK,+PKPo,)e J (^(;2_^2)[(;_^^)2_^2]- (9-10) 
From the above formulas we obtain 



1 



A{xi, X2){k2^ikiu - kl ■ k2gnv) 



^1^- - (47^)2 

, o/ .(PK-hPK-k2 . PK-h, PK-k2j \ 

+ B{xi,X2) I ^—^7 g^^u + PkijPku - _ h^PKu - J^T^^^^'^PkJ 

+ D{xi, X2) (^^i^^^T^d^i^ - ^^-^ki^hu + huPKv - , 

(9.11) 

where 



m 



K 



16a2e2{2[l - 2{xi + X2)]/i(^i^2) + xJi{z2) + X2[2/i(^|) - h{z2) + 1^)]} 
32a2e''{[2xl - xi{z + q)][-l2{zfz2) + h{zlz2)] 
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+ [2xiX2 -xi{z- q)/2 -X2{z + q) / 2][2h{zlzl) + 12(^1^2) - h{zlz2) 
- hizizl)] + [2x1 - X2{z - q)][l2{z^zl) - Hz^zl)]} 

4 / Tfl^ \ 

+ -a2e^ 1 + log ^ + (47r)2VMD^, (9.12) 




(9.13) 



D 2 

= 16026^/3 - 8026^/4 - 77026^(2; - q) log 

+ 16a2e2[2x2 - {z - q)/2][2h{ziZ2) - /i(^2)] 

+ 16026^(21/ - q)[h{zi) - A(l)/2] + 4a2e2[2xi - (^ + g)/2]/5 

+ (47r)2VMDi), (9.14) 

with the integrals given in the Appendix: 

4') = /' dz, dz2Z^z^' log (9.15) 

Jo Jo mt 



^^^^^ = / / d.2^, (9.16) 
mj(^ Jo Jo Di 



hm\^ I dzJ ''dz2Dilog^, (9.17) 
Jo Jo 

hm]c= I'dzMog^, (9.18) 
Jo 

h - /' dz^{Azl - Azi + 1) log (9.19) 
Jo mi 



and 
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D2 

Xi 



ml — 2ki ■ k2ZiZ2 — k\zi{l - zi), 

ml - klzi{l - zi), 
Pk -ki Pk- h 



m 



X2 



K 



m 



K 



(9.20) 



VMD^(a;i,a;2) = - 

V=0J,p 



Pk ■ (P 



K 



+ 



Pk ■ (Pk - ki) 



{Pk - ^2)^ - [pk - ki) 



m 



V 



(9.21) 



VMDBixi,X2) = - J2 Gvkrk2 

V=uj,p 



+ 



(Pk - - rri^ {Pk - h)^ - '^v 



, (9.22) 



VMBd{xi,X2)= J2 Gv 



ki ■ k2 



V=uj,p 



{Pk - kiY - m^' 



assuming the numerical values [54] 



(9.23) 



Gpm\ = 0.68 X 10-^ G^m\ = -0.28 x 10"^ 



(9.24) 



The loop calculation that we have just described provides all of the off-shell 
dependence scaled by the pion mass, and is of the form kf/ml. There can be an 
additional dependence of the form kf/A"^, where A fti 1 GeV. We cannot provide 
a model independent analysis of the latter. However, experience has shown that 
most of the higher order momentum dependence is well accounted for by vector 
meson exchange. Therefore we include the kl/A? dependence which is predicted 
by the diagrams of Fig. 28. One can recover the parametrization in ay neglecting 
the dependence on (p;^ — fci)^ and {pK~k2)^ in formulas 9.21-9.23, and performing 
the replacement [54] 
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Figure 30: Differential branching ratio dBR{KL 7r°7e+e-)/(i^ to 0{^^). 

where Ges ~ Gp+G^. This completes our treatment of the Kl — > 7r°7e"'"e~ amph- 
tude. 

The calculation we have presented in this chapter leads to the total branching 
ratio of 



BR{Kl TrV^e") = 2.4 x 10" 



(9.26) 



The decay distributions are presented in Figs. 30 and 31. 
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CHAPTER 10 



CONCLUSIONS 



Rare kaon decays are an important testing ground of the electroweak flavor 
tfieory. Witfi tfie improved experimental sensitivity expected in the near future, 
they can provide new signals of CP violation phenomena and an opportunity to 
explore physics beyond the Standard Model. 

The theoretical analysis of these decays is far from trivial due to the very low- 
mass of the hadrons involved. The delicate interplay between the flavor- changing 
dynamics and the confining QCD interaction makes very difficult to perform pre- 
cise dynamical predictions. Fortunately, the Goldstone nature of the pseudoscalar 
mesons implies strong constraints on their low-energy interactions, which can be 
analyzed with effective lagrangian methods. The ChPTh framework incorporates 
all the constraints implied by the chiral symmetry of the underlying lagrangian at 
the quark level, allowing for a clear distinction between genuine aspects of the Stan- 
dard Model and additional assumptions of variable credibility, usually related to 
the problem of long-distance dynamics. The low-energy amplitudes are calculable 
in ChPTh, except for some coupling constants that are not restricted by chiral sym- 
metry. These constants reflect our lack of understanding of the QCD confinement 
mechanism and must be determined experimentally for the time being. Further 
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progress in QCD can only improve our knowledge of these chiral constants, but it 
cannot modify the low-energy structure of the amplitudes. 

Addressing the two specific rare decays analyzed in this thesis, we can conclude 
that because of the three possible contributions to Kl i^^e^e" , the analysis of 
this process has multiple issues that theory must address. We have provided an 
updated analysis of all of its components. The goal of identifying direct CP viola- 
tion will not be easily accomplished. The decay rate by itself suffers from a severe 
uncertainty in the analysis of the mass matrix contributions. In the chiral analysis 
there is a free parameter, ^2, which is not fixed experimentally, and which has a 
strong influence on the decay rate. Measuring the related rate of Ks vr^e+e" 
would determine this parameter and will likely allow us to determine whether or 
not direct CP violation is present. 

Alternatively if the Ks branching ratio is not measured it may be possible to 
signal direct CP violation using the asymmetry in the e"*", e" energies. The direct 
and mass matrix CP violations have different phases, and the asymmetry is sensitive 
to their difference. There exist some combinations of the parameters where there 
remains some ambiguity, but for sizeable portions of the parameter space direct CP 
violation can be signaled by a simultaneous measurement of the decay rate and the 
energy asymmetry, as illustrated in Fig. 23. 

Progress in the theoretical analysis is possible, especially in the CP-conserving 
amplitude that proceeds through the two-photon intermediate state. Here both 
a better phenomenological understanding of the related decay Kl 7r°77, and a 
better theoretical treatment of the dispersive contribution should be possible in the 
near future. 

Overall, our reanalysis shows that the demands on the experimental exploration 
of this reaction are quite severe. The simple observation of a few events will not 
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be sufficient to indicate direct CP violation. The measurement of an electron 
asymmetry requires more events, and may or may not resolve the issue. Only 
the simultaneous measurement of Kg — > 7r°e'^e~ allows a convincing proof of the 
existence of direct CP violation. 

The behavior of the —>■ 7r^je'^e~ amplitude mirrors closely that of the process 
Kl — > 7r°77. The more complete calculation at order gives a rate which is more 
than twice as large as the one obtained at order E"^, despite the fact that the 
new parameter introduced at order E^ is quite reasonable in magnitude. This large 
change occurs partially because the order E^ calculation is purely a loop effect, while 
at order E^ we have tree level contributions, and loop contributions are generally 
smaller than tree effects at a given order. It was more surprising that the spectrum 
in Kl n^jj was not significantly modified by the order E^ contributions. These 
new effects are more visible in the low z region of the process Kl 7r°7e"'"e~. 

This reaction should be reasonably amenable to experimental investigation in 
the future. It is 3-4 orders of magnitude larger than the reaction Kl — > 7r°e''"e~ 
discussed above, which is one of the targets of experimental kaon decay programs, 
due to the connections of the latter reaction to CP studies. In fact, the radiative 
process Kl 7r^'je~^e~ will need to be studied carefully before the nonradiative 
reaction can be isolated. The regions of the distributions where the experiment 
misses the photon of the radiative process can potentially be confused with Kl — > 
7r°e"'"e~ if the resolution is not sufficiently precise. In addition, since the 7r° is 
detected through its decay to two photons, there is potential confusion related to 
misidentifying photons. The study of the reaction Kl — > 7r°7e+e~ will be a valuable 
preliminary to the ultimate CP tests. 



APPENDIX 



RELEVANT INTEGRALS 



In this Appendix we list the exphcit expressions for the integrals used in the 
calculation of chapter 9. We follow the notation of that chapter. For s < 4m^ and 
kf < Ami '^^ have: 



^ 4m^ — k\ arctan , 



\ Ami - kl 



s arctan , 



Ami ~ J ' 
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+ 



1 
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.2 
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^ Ami - kl 



(A.2) 



h{zl) 



kl 
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6^ 
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9 6a/^ 
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(A.3) 
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4 3m^ + kj 



^ (—Ami + 5mls — s^) 
+ ^ — , ^ ^ arctan , 



+ 



+ 



3s\/s\ Ami — s 



Ami — s 



mz 



2(s-A;?)2 



arctan , 



+ 2 



mi 



arctan 



4m^ 



ml 



'\ Ami ~ ^1 



+ 



ml 



kf 



3(s - A;f) 



2\2 



mt 



+ 



[—Ami + bmis — s 



s^/sx Ami - s 



arctan , 



Am^ 



{-Aml + bmlkl-k\ 



arctan 



kl^kl^Aml-kl 
mlkf 



kl 



{s - klf 



F{s)-F{ki) 



\ Ami - kj 

2mlkl 



[s-klY 



Ami — s 



arctan , 



4m?. — s 



^Aml - kl 



arctan 



kl 



\ Ami ~ kl 



(A.4) 



h{zl) 



+ 



+ 



+ 



2 



mz 



kl 



+ 



1 



9 A{s-kl) 2A{s-kl) 3(s-fcf)2 
ml {8ml — Qmls + s^) ' 
37 " 



ml - Smlkj - ^ 



arctan , 



6s\/s\/ Ami ~ ^ 

Sy/s^Aml 



Ami — s 



3(s - klf 



+ 



arctan , 



mti 



mz 



Ami 



kl kl^kl^Aml-kl 



arctan , 



mt. 



ml 



kl 



kl 



3{s-klf 

2 



s 



\ Ami ^ kl 

Iml — 6mls + s^) 

— arctan , 



m 



- + 



m 



2s^/s^jAml — s 
4 - 6mlkl + kf 



Ami 



kl 2kl^^Aml - kl 

m^kl 

{s — kiy ml ml 



■ arctan , 



kl 



\ Ami ~ kl 



s 2 J Ami - 

H ^1 arctan , 



Ami ~ ^ 



82 



+ — — arctan 



mi 



kl 



\ Ami ~ ^1 



-F{s)-rF{kl) 



mlkf 



F{kl) 



2^3 



F{s) 



3 J 4ml - kl 



arctan , 



arctan , 



4m?. — s 



kl 



^ Ami - kl 



(A.5) 



-^1(^1^2) 



+ 



+ 



+ 



13 ^ {-Qml 



^^)+i^!^l±^S^ arctan. 



144 24(s - kl) l2sJ~sJAml - s 



mz 



2{s - kl)^ 
2{s- klf 



F{kl) - F{s) 



ml (8m^ + 2mis - s^) 

+ ^ , - arctan , 



^•s &s^/s^ Ami ~ ^ 



ml (8m^ + 2mlkl - kf 



3^1 dkl^^^Aml - kl 
m^kl 

^J Ami ~ kl 



arctan , 



kl 



\ Ami - kl 



Ami ~ ^ 

7^ 



arctan , 



Ami — s 



arctan , 



kl 



\ Ami ~ kl 



Ami ~ ^ 



Ami — s 



(A.6) 



+ 



mz 



2{s - kf) (s - kl)^ 
kj \ \/Aml - s 



F{s) - F{kl) 



(s-kl)^ 
^jAml - kl 



arctan , 



4m?. — s 



arctan , 



kl 



^ 4m2 — kl 



(A.7) 



^2(^1^!) 
m^ 



'kl 



6{s-ki) {s-kiy 



m: 



kl)' 



2^ Ami ~ arctan , 



4m?. 



83 



- 2 J Ami - ki\ ki arctan . - — „ ^ , „ + 
^ ^ \ 4m^ — kf 



+ 



+ 



2kf 



3{s-kiy 



ml (—Ami + bmls — s^) 

— + ^ — - arctan . 

"S sJs\ Ami - s 



Am^ — s 



ml {-Am^ + hm^kl - kf] 



arctan , 



ki kl^^^Aml - kl \ - kl _ 



kf 



2mlkl 
{s - klY 



F{s) - F{kl) - 



2. Ami - s 



arctan , 



Ami — s 



2 J Ami - kl 



arctan 



kl 



^ Ami ~ kl 



(A.8) 



mji 



+ 



+ 



+ 



1 

~12{s-kl 

mt 



Ami 



_1 f 6ml + kl 

- klf \ A 
I2mlkl-k\\ 



mi 



{s - klY 

k^^^i,nl-k 



A;? s^/sJAml-s 

— H — : arctan 



mi 



m„ 



Ami — s 



mt 



kl 



{s-kl 



2\A 



m; 
s 



■ arctan 



5m: 



kl 



\ Ami - kl 4.m^ 

Gmls + s^) 

arctan , 



J4_ 

Ami 



s^/sJAml - s 



Ami 



mi 



Qmlkl + ki] 



kl kl^^^Aml - kl 
2^ Ami - ^Vs 



arctan , 



k^ 



\ Ami ~ kl 



Smjkl 

{s - klY 



mi 



arctan 



Ami ~ ^ 



mi 



kl 2^Aml-kl^kl 



mi 



mi 



arctan 



kl 



\ Ami ~ kl 



- F{s) + F{kl) 



3mlkf 

is - klY 



3 J Ami ~ s 
F(s) p arctan 



Ami ~ s 



- F{kl) + 



3 J Ami kl 



arctan , 



kl 



\ Ami ~ k\ 
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+ 



Q{s-kl) {s-klY 



Ami 



kf^Aml - kl 



arctan , 



kf 



^ Ami - kj 



77fv4<-^? arctan ^^-^ 



{2kjml + kfs — Qmls) 



arctan 



4m?. — s 



(A. 10) 



24(s - fcf) 12(s - kf)^ 
kf 



{6ml - kl) - 



ml 



[s - klY 



F{s) - F{kl) 



{s - kl)^ 



(8m^ + 2m^s — s^) 
6sy/s^Aml — s 



+ 



arctan 



Ami ~ ^ 



ml {Smj + 2mlkl - kf] 

3^1 6kl\/kf^/Aml~kf 

2mlkl 



arctan 



\ Ami - kl 



kir 



Ami ~ kl 



arctan 



kl 



\ Ami ~ kl 



■ arctan 



Ami ~ ^ 



(A.ll) 



+ 



12{s-kl) 2{s-kiy 



Ami 
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— {2khtr + /,-r,s — (y/zr s) 

4m^ — s '■ — ^ ^ '■ — arctan , 



Ami — s 



(A.12) 



+ 



+ 



+ 



13 



13 



-12^^^+144 



ml 



kf 



2{s - 



(8m4 + 2mlkl - kf) 



F{s) - F{kl) 



arctan , 



QkiJkiJAml - kl 



kl 



\ Ami — kl 



(8m^ + 2mls — s^) 
Qsy/s^Aml — s 

klml 



arctan , 



4m^ 



s-kl 



Ami. ~ ^1 



arctan , 



kl 



^ Ami - ^1 



Ami 



arctan , 



^/s Y Ami ~ =^ 

-{2klml + kls - lOmls + s^) 



12s 



arctan , 



Ami 



(A.13) 



hm^ = — 



5kl Ami {^rnl-kl)^Aml-kl 



+ 



arctan 



^ Ami ~ ^1 ' 



(A. 14) 



8 Sml 2{Aml-kl)JAml-kl 

-- + -:7TY — arctan 

9 3ki 



kl 



\ Ami ~ ^1 ' 



(A.15) 



In the cases when s > Ami > 4m^, we have to perform the substitutions: 



arctan 



4m?. — s 



arctan , 



kl 



Ami - ^1 



1 

Ti 
1 



log 
log 



Ami/ s 



l + ^l-Aml/s^ 



l-Jl- Aml/kl 



+ m 



+ m 



1 + Jl- Aml/kl 



(A.16) 
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and 



Ami — s ~^ iys — 4ml, 
^jAml - kl i^kl - Ami, (A.17) 



respectively, in formulas A.1-A.15. 
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